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PLOT 


Calculation of Linear Regression Line, Plotting of Data Points 
and Plotting of Linear Regression Line 


DECUS Program Library Write-up DECUS No. 8-202 

SCALE DETERMINATION 


The computer will calculate a scale for the x axis and the y axis from the minimum and 
maximum x and y values, respectively. Therefore, one must enter as the first two "data 
points" value for x and y that will determine the x and y scales. These two points will 
be pidtted on the graph paper, but will not be included in the calculation of the linear 
regression line. 


Example One 

If you want to plot on a half page (8 1/2 x 11) x values which are between 1.0 and 8.0, and 
y values between 1.0 and 5.0, then it would be convenient to add the two points 0,0 and 
9,6. The X scale would then be typed out as 0.1 and the y scale 0.1, indicating that each 
little square on the graph paper corresponds to 0.1 unit. The x and y axis values for the 
heavy graph lines would be 1.00, 2.00, 3.00, 4.00, 5.00, etc. 

Example Two 


Say the x values on the half page are between 3.2 and 4.5. Good beginning and end points 
would be 3 and 4.8 respectively. The ten lines on the x axis should be marked: 3.00, 
3.20, 3.40, 3.60, 3.80, 4.00, 4.20, 4.40, 4.60, 4.80. 


General rule: find two end points, between which the actual data points will fall, and so 
that the difference between the two end points can be divided by the number of large squares 
(listed below) on the corresponding axis giving a convenient scale on that axis. 


Half page x axis contains 9 large squares 
Half page y axis contains 6 large squares 
Full page x axis contains 14 large squares 
Full page y axis contains 9 large squares 

The actual markings of the axes can be done with the Alphanumeric Plotter Program, ALPHA, 
DEQJS No. 8-203, or automatically when an 8K memory is available. , 

LOCATIONS 

The program uses the following locations: 

5, 7, 40-62 for the floating point package 

2, 3, 15, 16, 17 for all sections of the program 

20-177 for all sections of the program 

200-2453 , 4650-4751 for all sections of the program 

2460-3553 for storage of x values, 3 locations per data point 

3554-4647 for storage of y values, 3 locations per data point 









A total of 190 x-y combinations can be stored. 
OPERATING INSTRUCTIONS 


1. Load binary tapes for this program in the usual manner. Also, the Floating Point Package, 
No. 3 for EAE Type 182, Digital 8-25-U BIN (Pels Tape ^ 7 ) should be loaded. When ex¬ 
tended memory is available, the Alphanumeric Plotter Program should be loaded in the basic 
memory (Field 0) and PLOT and the Floating Point Package both in Field 1. After the points 
and linear regression line are plotted, the computer will then go to the Alphanumeric Program 
automatically. 

2. Set RSW to 0200 in Field 1 (or Field 0, if no extended memory is available). "LOAD 
ADD" and "START". 


3. For a half page, put bit 10 up. For a full page, put bit 10 down, 
are the same as used in the Alphanumeric Plotter Program. 


These switch settings 


4. When a print-out of the entered data is required, put bit 1 up. When no print-out is needed 
put bit 1 down. When entering the data by hand from the keyboard, choose the print-out 
option. When entering the data from a paper tape, not typing out will speed up the operation. 


5. When the user does not want the regression line to be drawn on the graph, put bit 9 up. 
The computer will then calculate the regression line, type out the values, type the x scale and 
y scale, and plot the points. Then it will go to a next regression line calculation. Put bit ^ 
up when automatic scaling is required. 


6. Before entering the data, set the plotter pen manually exactly on point A, indicated in Figure 1 
or Figure 2, depending on whether the user wants a half page or a full page format. Location 1320 
contains 5331 which causes the computer to bypass a section of the program which automatically sets 
the pen on point A. We found that the pen might be off-center and so does not go exactly to point A. 
Therefore, we inserted the JMP. The user is invited to try the subroutine by changing 1320 to 7300^^ 
It is possible that the counters MIOOO and N18 need adjusting. 


7. To enter the data from the keyboard or to punch a paper tape; 

"RETURN" or letter y (see explanation at the end of this paragraph). 

X SPACE y SPACE These values must be the additional points used 

o o 

x^ SPACE y^ SPACE to calculate the scales, as explained above. 

X 2 SPACE y 2 SPACE This is the first actual data point. 

X SPACE y SPACE This last y-value must be followed by RETURN (EOL). 

When the x values are to be plotted on the x axis and the y values on the y axis strike or punch 
"RETURN". When the y values are to be plotted on the x axis and the x values on the y axis, 
strike or punch Y. 
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8. When the return is sensed by the computer it will calculate from the data points the linear 
regression line and type out the values: 



In the formula y = Bx + C ± standard error 


SE J 


See last page for the 

PCT 

Prediction error in percent 

formulas used in the 

RSQ 

?■ 

calculations. 

R 

r (correlation coefficient) 


SER 

Standard error of r 


N 

Number of entries (x-y combinations) 



9. The computer then waits. Now any value of x can be entered followed by a terminator (see 
below) and the computer will calculate and type out the corresponding y value and the y value ± 
standard error. Selecting x values close to either end of the graph will allow the user to draw the 
band around the regression line very easily. For this, two x values will be required; terminate the 
first X by a space. This will permit the computer to accept another x value. The last x must be 
terminated by either a RETURN or a decimal point. 

RETURN Will bring the computer back for another regression line calculation. 

DECIMAL POINT Will start the plotting of the points and the regression line plot. 

SPACE Will bring the computer back for another x value. 

10. If the X value is an integer, a following decimal point will be interpreted as a real decimal 
point. In order to terminate the entry, a second decimal point is required. Before striking the 
terminating decimal point, establish the proper switch setting: 

a. Bit 0 down: The plotter pen will stop after the regression line is drawn, 

ready to start in the Alphanumeric program. 

b. Bit 0 up: The plotter pen will go all the way to the right on the x axis 

after the present plot is completed, and start drawing the x and 
y axis lines and the axis notations with the Alphanumeric Plotter 
Program (see 1). 

11. When the decimal point is used as terminator, the computer will first calculate the y and y ± 
standard error, rank the x values in descending order and then calculate the x and y scales and type 
these out. Then it will start plotting, first two end points and then the actual data points. 

12. Figure 1 is a completed plot in half page format. Figure 2 gives the same points in full page 
format. Figure 3 is a typical typ>e-out (for the half page option of Figure 1). 

13. Resume of switch settings: 

Bit Up Down 

0 After plotting the points and regression After plot is finished, prepare 

line, plot the scales automatically. for Alphanumeric plot. 






Down 


Bit Up 

1 Print-out of entered data 

9 Do not plot regression line. 

10 Half page format. 


No print-out of entered data. 
Plot regression line. 

Full page format. 


14. When only the basic memory is available, the two programs should be used one after 
the other. In this case, location 1721 in this program could be changed to HIT (7402). 

This will halt the program. 

15. In some cases, a slight job will appear at certain intervals in the regression line. This is 
due to a limitation imposed by the size of the plotter motor step and is more noticeable with the 
felt pen than with the ball point pen. 
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Regression and Correlation 


y = bx -j- c+ standard error 


b = NZxy - (Zy) (Zv) 

NZy^ - (Ex)2 



bUx 

N 


standard error = 


y Z(y calculated - y experimental) 

'' N 


2 


= [NZyy - (Zx) (^v)]“ . 

tNZx^ - (Zx)2j InZ y2 - (Zy) 


r = Vr‘ 


standard error of r = 1 - r' 


VN 


PCT = 


\ /Z /y calculated - y experimentalV 
\ / \ y calculated y 

V N-2 
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Figure 3. Type-out for the four sets of data on the 
test tape for half page format. 

B +)2f.6666666E+J3rj2f 
C +0r.476 8371E-J2f6 
SE +jar.4396219E-0r6 
PCT +0r.l762792E-j3r4 
RSQ +0.100g(0gl0E+gll 
R +9f,19!09!9l9l0TS+01 
SER +gf.000gf000E+i2f0 
N +0.5000000E+01 

+0f.6666671E+jafjar +0.6666666E+00 +1?. 6666674E+!3^0^ 

X SCALE +jar.9999999E-0fl 
Y SCALE +0.9999999E-0fl 


B +j3r.9 71427 4E+J2f)2f 
C -!2f.lj2r2 8566E+J2fl 
SE +(?. i385iar5j2rE+!2r(2r 
PCT +!ar. 234(2r5!3r8E+!ar2 
RSQ +j2f.99 49 423E+|3r!2f 
R +0.9974679E+|y(3f 
SER +12r.l91162 8E-0r2 
N +0.7000000E+01 

+|2f.3342856E+0fl +)2f. 32!af435|2fE+jafl +!3r. 348136jafE+(3fl 

X SCALE +jar.9999999E-!3ri 
Y SCALE +j3r.9999999E-|3fl 


B +0 . 1000000E+01 
C -0^1999999E+01 
SE +!ar.3693562E-|ar6 
PCT +J3r.288675j2rE+!3r2 
RSQ +0.0000000E+00 
R +0.1000000E+01 
SER +0.0000000E+00 
N +0.5000000E+01 

+!y.25j2rjaf!3riafj3fE+jari +!af.25|2r0f)2fj3r0rE+0ri +j3f.25j2f|2ri2rjar0fE+)2ri 
X SCALE +)ar.9999999E-!3ri 
Y SCALE +jy.9999999E-01 


B +0.7500000E+00 
C +0.5000000E+00 
SE +0.0000000E+00 
PCT +0.0000000E+00 
RSQ +0.1000000E+01 
R +0.1000000E+01 
SER +0.0000000E+00 
N +0.4000000E+01 

+|2f.425|2fWE+0^1 +)2f.425jar0fJ2fi2fE+(2ri tj?. 425J2fWE+!2ri 
X SCALE +1^.999999 9 E-jafl 
Y SCALE +j 3 r. 9999999 E-j 3 ri 
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